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The observation of fluctuations of parity-odd angular observables at RHIC has been interpreted 
as a signal of a local parity violation. We show how the observed correlations can be explained by 
local charge conservation at freeze-out combined with elliptic flow. Calculations from a blast wave 
model, which overlays thermal emission onto a collective flow profile, are shown to account for the 
experimentally observed signal. 
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The prospect of observing parity violation from the 
strong interaction in relativistic heavy ion collisions [1 
has recently gained great attention. Measurements of 
the STAR collaboration have revealed large fluctuations 
of parity-odd observables and the signal was proposed to 
originate from a local parity violation in QCD leading to 
opposite directions of preferred emission for oppositely 
charged particles [2]. Since the direction of preferred 
emission fluctuates randomly, the studies are confined 
to charge-dependent azimuthal angle correlations. The 
observable measured by STAR is 
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where a and /3 represent positive or negative charge. 
Ma and Mp are the corresponding multiplicities, the az- 
imuthal angles (j) are measured about the beam axis, and 
the reaction plane is at (/) = without loss of general- 
ity. In the case where a and /3 refer to the same charge, 
the i = j terms are excluded from the sum. The pro- 
posed source of the parity-odd fluctuations [T] is an elec- 
tric field generated perpendicular to the reaction plane 
of the initial collision. Whereas the passing ions would 
be expected to create a non-zero coherent magnetic field 
for events with non-zero impact parameter, there should 
per se be no coherent electric field. However, through the 
anomalous coupling to the parity-odd topological charge 
in the QCD sector, Ea • Bq, a coherent electric field can 
be generated, whose direction would fluctuate from being 
parallel to being anti-parallel with the coherent magnetic 
field. For each individual nucleon-nucleon collision the 
sign of E(j • Ba within each individual flux tube would 
vary randomly, thus making the generated electric field 
coherent only in as much as there are several times more 
particles than fiux tubes. The phenomenon is referred to 
as the "chiral magnetic effect". Since one is looking for 
effects from a coherent electric field that is out-of-plane, 
the angular correlation should be constructed so that it 
compares correlations between same-sign and opposite- 
sign pairs, and so that it compares out-of-plane vs. in- 



plane correlations. The observable of interest is then 
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where it has been assumed that there are equal numbers 
of positive and negative charges. This is approximately 
true for the energies at RHIC. Then M = M+ + M_ is 
the total charged multiplicity, and "ss" refers to the cor- 
responding sum with same-sign pairs. A coherent elec- 
tric field out-of-plane (0 = 7r/2) would give a positive 
correlation between the sin (f> terms for same-sign pairs 
and a negative correlations between the sin terms for 
opposite-sign pairs. Thus, 7p would be positive, and is 
indeed observed to be positive. By subtracting the in- 
plane correlations, i.e., the cos (j)i cos 4>j terms, the intent 
is to eliminate many other sources of angular correla- 
tions that are not dependent on the average direction of 
the pair. Although the chiral magnetic effect would ex- 
plain the sign of the data, the magnitude of the observed 
signal is several orders of magnitude higher than some 
estimates Oil]. 

However it has been pointed out that local charge con- 
servation at breakup combined with elliptic flow might 
explain a major part of the observed signal [31. Simi- 
larly cluster particle correlations have been proposed as 
another alternative explanation [5] Charge conservation 
gives positive correlations between opposite-sign parti- 
cles resulting in positive terms for cos (pi cos (pj as well 
as sin (pi sin (pj (opposite to the parity signal) . But since 
there are more such pairs in-plane than out-of-plane due 
to elliptic fiow, and additionally because their correla- 
tions are stronger in-plane than out-of-plane (also oppo- 
site to the parity signal) the expression in ^ is dom- 
inated by the cos (p terms for opposite-sign pairs. Thus 
this effect also results in positive values of 7p. In this let- 
ter we investigate how such correlations can reproduce 
the observed signal by considering a simple blast wave 
model in which particles are generated according to ther- 
mal sources moving with collective flow. By constraining 



each thermal source to emit equal numbers of positive 
and negative particles, one can generate non-zero values 
of 7p. The goal of this letter is to quantitatively ex- 
plore the degree to which the measurement of 7p can be 
reproduced with a simple thermal model incorporating 
collective flow and local charge conservation. 

In all fundamental processes charge is created in bal- 
ancing pairs that are produced at the same point in space 
time. When the motion is highly collective the correla- 
tion in space time translates to a tight correlation of bal- 
ancing charges in momentum space. This correlation is 
strongest when charge production takes place late in the 
collision or diffusion is small [6] . A differential observable 
that has been exploited for measuring such correlations 
is the charge balance function [5], which measures the 
chance that a charge at angle </> has a balancing charge 
emitted with angle </> -I- Acf). Results from the STAR col- 
laboration have shown that balancing charges are likely 
to be emitted in a narrow range in rapidity [3| , consistent 
with balancing charges being emitted from small neigh- 
borhoods such that the relative rapidities are determined 
mainly by the thermal motion at breakup [5]. Similar 
results have recently been generated for azimuthal an- 
gles [9j. Additional contributions to the charge balance 
funtion due to small-angle correlations from final-state 
interactions have been considered in j8j, but only affect 
the results for jp at the level of a few percent and will 
therefore be negelected. For equal numbers of positives 
and negatives, the charge balance function is defined by 
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where Na^p{(j), Ac/)) is the number of pairs where the type 
a was emitted at (f) and the type /3 was emitted at (jj+Acj). 
The like-sign subtraction statistically isolates the balanc- 
ing partner. The same analysis could be done in the con- 
text of correlations, which for this case are identical to 
balance functions aside from an extra factor of dM/d(f>- 
The correlation 7p can be expressed in terms of mo- 
ments of the balance function, after being combined with 
the angular distribution dM/d(f>, 
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where one has inserted the definition of the balance func- 
tion ([3]) into the definition for the parity observable ([2]), 
and used the angle addition formula, cos(0i + 0j) — 
cos(20i) cos(A0) — sin{2(j)i) sin(A0). 

To understand the degree to which charge conserva- 
tion affects 7p, a thermal blast wave model was modi- 
fied to incorporate local charge conservation. Blast-wave 



models are simple parameterizations of the breakup con- 
figuration. For this study, the model used by STAR to 
fit elliptic flow data and spectra was employed p5] . The 
model parameters are the breakup temperature Tkin; the 
maximum collective velocities in the in-plane and out-of- 
plane directions, and the spatial anisotropy of the elliptic 
shape. The elliptic anisotropics were chosen to flt the el- 
liptic flow observable V2, which quantifies the degree to 
which more particles are emitted in-plane than out-of- 
plane. 
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The anisotropy is driven by the elliptic shape of the initial 
fireball, as viewed transverse to the beam. The quantity 
V2 is a staple of RHIC science and has been analyzed both 
as a function of the centrality of the collision (for very 
central collisions the initial shape has little anisotropy), 
as a function of transverse momentum (the anisotropy 
is stronger for higher pt particles) and as a function of 
species type (more massive particles are more sensitive to 
collective flow) [TU]. By fitting both spectra and V2, blast 
wave parameters were determined for several centralities 
by STAR in [TU]. In this parametrization the breakup 
configuration is characterized by four parameters: the 
kinetic freeze out temperature, the transverse rapidties 
in-plane and out-of-plane and the spatial anisotropy. The 
usual method to apply a blast wave model would be to 
choose a collective velocity consistent with the blast wave 
parameterization, generate a particle according to the 
thermal distribution characterized by the temperature 
and collective velocity, then repeat for several particles. 
For this calculation, an array of particles is generated 
rather than a single particle. The array is chosen consis- 
tent with a canonical ensemble with a fixed volume (64 
fm'^) and a chemical temperature, Tchem = 175 MeV [TT] . 
so that electric charge, strangeness and baryon number 
all sum to zero. The particles are then individually as- 
signed momenta according to the (kinetic) breakup tem- 
perature and collective flow. In this way, charge conserva- 
tion is enforced in the most stringent way, with balancing 
charges being emitted from the same source velocity. If 
the charges had been created early and diffused before 
most of the collective flow developed, balancing particles 
might have been emitted from regions with significantly 
different collective velocities. Analysis of charge balance 
functions as a function of (mid-)rapidity suggest that the 
emission of balancing charges is indeed highly localized 
for central events [S], but less so for peripheral events. 
Given that the collective velocity gradients in the trans- 
verse direction are smaller, the effects of diffusion are 
expected to be smaller for balance functions binned as 
a function of azimuthal angle, which makes constraining 
the balancing charges to originate from sources with the 
same collective velocity reasonable. 

Balance functions B{(j), Ap) from the blast wave model 



described above are presented as a function of At/) in Fig. 
[ijfor events with centralities of 40-50% using STAR's pa- 
rameters [10:. Here 0% centrality corresponds to zero 
impact parameter and we refer to [lOj for more details 
on the classification. The balance function for = 0° 
(in-plane) is narrower than the balance function for (p = 
90° (out-of-plane) . The stronger focussing of balancing 
charges derives from the greater collective flow in-plane 
vs. out-of-plane. For = 45°, the distribution is biased 
towards negative values of Atji. This is expected given the 
elliptic asymmetry, V2 > 0, which leads to more balanc- 
ing particles toward the 4> = 0° direction as opposed to 
4> = 90° . Depending on which quadrant cj) is located, the 
balancing charge tends to be found more towards = 0° 
or (p — 180°. The lower panel shows the moments of 
5(0, A0), 
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where 
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is the normalization of the balance function and repre- 
sents the probability of detecting the balancing charge 
given the observation of a charge at 0. It would be unity 
for a perfect detector, but is reduced by both the finite 
acceptance and efficiency of the experiment. The quan- 
tity Cb{4>) determines the width of the balance function 
and would be unity for a very narrow balance function 
whereas it vanishs in the case where the balancing charges 
were emitted randomly. The quantity Sf,(0) measures the 
degree to which the balance function is asymmetric under 
a reflection symmetry of A0 — > — A0. For pairs around 
= 45° this corresponds to the probability for the bal- 
ancing charge to be emitted in in-plane direction vs. in 
out-of-plane direction. 

The parity observable jp can be expressed, using the 
moments of the balance function defined in m , 
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where we introduced 



V2c 
V2s 



{fm ^ 



(Cfc(0)cos(20)) - W2(C6(0)), 

(sb(0)sin(20)), 
^/d0^..(0)/(0). 



(9) 



The three contributions to 7p derive from: a) having 
more balancing pairs in-plane than out-of-plane (w2(cfc)), 
b) having the in-plane pairs being more tightly correlated 
in A0 than the out-of-plane pairs (w2c)and c) having the 
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FIG. 1: (color online) UPPER PANEL: Balance function for 
40-50% centrality shown as function of the relative angle in- 
cluded by balancing partners for = 0° (black squares), 45° 
(red triangles) and 90° (blue circles). LOWER PANEL: The 
moments of the balance function, Ci,(0) and S6(0), represent 
averages of cos(A0) and sin(A0) across the balance function. 
These are plotted as a function of for various centralities. 
The structure of C(,(0), which is maximized at = 0°, 180°, 
illustrates how the balance function is narrower for in-plane 
emission and more central collisions, while the structure of 
Sbifj}), which is positive for = 135°, 315° and negative for 
= 45°, 225°, shows how balancing charges prefer to be emit- 
ted in the in-plane direction. The oscillations increase for 
more peripheral collisions. 



balancing charge more likely being emitted towards the 
event plane {v2s)- The first term was estimated in |3] to 
be a significant fraction of the observed signal, but the 
latter two terms could not be estimated without a more 
detailed model like this one. The contributions for 7p 
obtained from the blast wave calculation are displayed in 
Fig. [2] For better visibility the values are scaled by the 
multiplicity (for the STAR data we use the experimental 
multiplicity 12J to account for efficiency and acceptance). 
To compare to the STAR data [2] we use the same accep- 
tance cuts in transverse momentum and pseudorapidity 
and we assume perfect detector efficiency for the blast 
wave calculation. The necessary efficiency correction is 
done by rescaling the results to reproduce the experi- 



mental normalization of the balance function ^, i.e. we 
multiply the expressions for {cb{4')), vic and V2s by the 
ratio of experimental to blast wave normalization. 
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FIG. 2: (color online) Parity observable from STAR (black 
stars) and blast wave calculations (red squares). The three 
contributions to the signal are defined in Eq. ([8| and are 
plotted with dashed lines. V2{cb) (green circles) derives from 
having more balancing pairs in-plane than out-of-plane while 
V2,c (blue triangles) quantifies the degree to which in-plane 
pairs are more tightly correlated than out-of-plane pairs. «2,s 
(orange triangles) reflects that the balancing charge is more 
likely to be found towards the event plane. 

For central collisions the STAR data is well reproduced 
by the blast wave model, however for more peripheral 
collisions the model produces higher correlations than 
experimentally observed. As we assumed that balanc- 
ing particles are emitted with the same collective flow, 
the correlations from the blast wave model shown in Fig. 
[2J are the strongest possible contribution to the signal in- 
duced by charge conservation. Consequently the model 
is expected to over-predict the measured signal. Whereas 
for central collisions this assumption has been shown to 
be consistent with measured balance functions [5], the 
locality of charge conservation at breakup appears to be 
less exact for peripheral collisions. The STAR analy- 
sis of charge balance functions [5], shows a significant 
broadening of charge balance functions for larger impact 
parameters, a part of which is expected from the higher 
kinetic freeze-out temperature and less transverse collec- 
tive flow. By relaxing the stringent conditions of having 
the balancing charges emitted from the same azimuthal 
angle, the balance functions will become wider and the 
accompanying correlation jp should move down toward 
the data. In the limit that the balancing charge were 
allowed to come from any point in the blast wave, the 
contribution to jp would vanish. 

The calculations presented here demonstrate that local 



charge conservation overlaid with elliptic flow can read- 
ily explain the difference between the opposite-sign and 
same-sign correlations seen by STAR. The over-simplified 
picture of the decoupling provided by the blast wave 
model is rather crude and in some ways the parame- 
terization is arbitrary, but should be accurate to within 
a few tens of percent. To better model how the ob- 
served 7p is induced by these mechanisms, one should 
employ a more detailed model of both the collision dy- 
namics and decoupling and of the correlation in space- 
time between balancing charges. We emphasize that 
the mechanism described here only explains the differ- 
ence between the opposite-sign and same-sign correla- 
tions. This is sufficient to show that large parity fluctua- 
tions are not warranted by the data, although the strong 
same sign-correlations observed at STAR [2] have yet to 
be explained. It was suggested in [4^ t^i^t these might 
be induced by momentum conservation, however there 
is so far no quantitative theoretical estimate of this ef- 
fect. A more completely satisfying model would describe 
both the same-sign and opposite-sign correlations inde- 
pendently, and would also reproduce the more detailed 
differential correlations, such as charge balance functions, 
in addition to the integrated correlations, jai3- We ex- 
pect such analyses to be pursued in the near future. 
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